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Abstract 

Scaling similarity solutions of three integrable PDEs, namely the 
Sawada-Kotera, fifth order KdV and Kaup-Kupershmidt equations, 
are considered. It is shown that the resulting ODEs may be written 
as non-autonomous Hamiltonian equations, which are time-dependent 
generalizations of the well-known integrable Henon-Heiles systems. 
The (time-dependent) Hamiltonians are given by logarithmic deriva- 
tives of the tau- functions (inherited from the original PDEs). The 
ODEs for the similarity solutions also have inherited Backlund trans- 
formations, which may be used to generate sequences of rational so- 
lutions as well as other special solutions related to the first Painleve 
transcendent. 



1 Introduction 



The six Painleve transcendents (PI- VI) have had a considerable amount of 
attention devoted to them in recent years, and have been studied from many 
different points of view. Their original discovery came about from Painleve's 
classification of second order ODEs (of a particular form) having no movable 
critical points. They have also been approached by way of isomonodromic 
deformation of linear differential equations |L4j] , or via abelian integrals and 
algebraic geometry [ |25[| . Furthermore, they have found numerous physical 
applications. For example, the first Painleve transcendent (PI), 

d 2 w 9 

^ = &»»+,, (i.i) 

arises as a string equation in the matrix models of 2-D gravity (T^j, while 



other Painleve transcendents govern the behaviour of correlation functions 
in certain models of statistical mechanics ||, quantum field theory |2!| and 
random matrices [0 (see chapter 7 of |3J for a more complete survey). Higher 
order ODEs of Painleve type also occur in many situations, but as yet there 
is no complete classification of such equations. 

Another context in which ODEs of Painleve type arise naturally is as 
similarity reductions of integrable PDEs. A great deal of the renewed interest 
in Painleve equations was generated by the conjecture of Ablowitz, Ramani 
and Segur Q, which states that all similarity reductions of integrable PDEs 
are of Painleve type. For example, the modified Korteweg-deVries (mKdV) 
equation, 

v t = v xxx - 6v 2 v x , (1.2) 
has scaling similarity solutions of the form 

v(x, t) = (-3t)~*y(z(x, t)), z= (-3t)~sx, 

where y(z) satisfies 

y" = 2y 3 + zy + a, (1.3) 

for some constant a (' denotes -p), which is the second Painleve equation 
(PII). An important feature of the similarity reductions of soliton-type equa- 
tions is that they inherit much of the stucture associated with integrability, 
such as Backlund transformations and solutions in terms of special functions. 
In the work of Okamoto P71, each of the equations PI- VI is written as a 



non-autonomous Hamiltonian system, with a polynomial Hamiltonian, e.g. 
(|1.3[ ) may be obtained from the Hamiltonian 

h = \p 2 -{y 2 + l)p-{a + \)y (1.4) 
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(where y, p are the canonically conjugate variables, and z denotes time). 
Each of the Hamiltonians used in p7| has explicit dependence on the time 
z, and is obtained from a holomorphic tau-function, 

h(z) = -^-log[r(z)]. 

The tau-function provides a simple way to encode both the solutions of the 
equations and their Backhand transformations. In the case of PII, this tau- 
function essentially coincides with the tau-function of mKdV/KdV (up to 
a factor of exp[^ 3 /24]). Manin |25| has found an alternative Hamiltonian 
description for PVI, involving the Weierstrass p function, but his approach 
is difficult to generalize to higher order equations of Painleve type. 

In what follows we apply the approach of Okamoto, using polynomial 
Hamiltonians and tau-f unctions, in order to analyze the the scaling similarity 
solutions of three different integrable, fifth order evolution equations, namely 
the Sawada-Kotera (SK) equation, 

ut = u 5x + 5uu 3x + •")//,//,,,. + 5u 2 u x , (1.5) 

fifth order KdV, 

fh = u 5x + 10uu 3x + 20u x u xx + 30u 2 u x , (1.6) 

and the Kaup-Kupershmidt (KK) equation, 

u t = u 5x + 10uu 3x + 25u x u xx + 20u 2 u x . (1.7) 

The ODEs for scaling similarity solutions of each of the flows in the SK, 
KdV and KK hierarchies, as well as their modified hierarchies (related by 
a Miura map), are presented in section 2. In the third section, by a di- 
rect extension of Fordy's results [T!| on stationary flows, it is shown that 



the scaling similarity solutions of the fifth order equations in each hierarchy 
are related to non-autonomous generalizations of the integrable Henon-Heiles 
systems. We describe the tau-functions, Miura maps and Backhand trans- 
formations associated to these systems, and use this machinery to generate 
sequences of special solutions in section 4. The essence of our method is 
that the Hamiltonian formalism for stationary flows [|17]] may be extended 
to other (non-stationary) similarity reductions. We briefly discuss further 
developments in our concluding section. 

Note. In their most basic form, these results were originally presented 
in 1 19], although these proceedings do not seem to be generally available 
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at present; my thesis |2(J contains a more thorough discussion. Shortly 
after my thesis was submitted, I was made aware of Kudryashov's paper f24f 
concerning similarity reductions in the KdV/mKdV hierarchy. The results in 



24] have some overlap with my own, but do not make use of the Hamiltonian 



formalism. 



2 Scaling Similarity Solutions in the Sawada- 
Kotera, KdV and Kaup-Kupershmidt Hi- 
erarchies 

In this section we develop a concise notation to describe the scaling similarity 
solutions of any one of the flows of the SK, KdV or KK hierarchies. We then 
apply it to some particular examples, including the similarity solutions of the 
fifth order equations which we study in detail in sections 3 and 4. 

2.1 General Description of Scaling Similarity Solutions 

SK and KK have only one Hamiltonian structure, but KdV is bi-Hamiltonian, 
and here we will be using the second Hamiltonian structure. Following Fordy 



16 , we are able to consider all three hierarchies at once. The nth flow in 



each of the hierarchies can be written as 
du 

— = (8l + 8aud x + Aau x )5 u H n [u}, (2.1) 
ot n 

where a = 1/2 for SK and KdV, a = 1/4 for KK, and H n is the nth Hamil- 
tonian density for the hierarchy in question. There is also a Miura map from 
the modified versions of the hierarchies, given by 

u = -v x - 2av 2 =: M[v]. 

The nth flow of the modified hierarchy is 

= {-d x )5 v H n [M[v}}. (2.2) 

ut n 



The Miura map means that given v satisfying ( |2.2|) for each n, the corre- 
sponding u = —v x — 2av 2 satisfies (|2.1|). Further details on these hierarchies 



may be found in [15 



The nth flow of the hierarchy is unchanged by the scaling 

x — > (3x, t n — > p m t n , u — > j3~ 2 u, 
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where m = m(n) is a scale weight dependent on the hierarchy and on which 
flow is being considered. Similarly the modified flow is invariant under the 
same scaling but with 

v (3~ l v. 

Hence there are scaling similarity solutions of the form 

u(x,t n ) = 6 2 (t n )w(z), z = x6{t n ), ^- = 9 m+1 . 

dt n 



(We use the notation of |IT|.) The corresponding similarity solution for the 
modified flow is 

v = 0(t n )y(z), 
with the scaled Miura map giving 

w = -y' - lay 1 =: M[y] 

(' denotes ^ throughout). 

In the context of an integrable hierarchy, it is customary to think of the 
dependent variable as a function of all the times, 

U = = X,t 2 ,t 3 , ...). 

For the scaling similarity solutions of the nth flow, it is necessary to drop 
dependence on anything other than x and t n , and it is inconsistent to consider 
the other flows simultaneously. To incorporate the other flows extra similarity 
variables must be included (e.g. in similarity reductions of the mKP equation 
[|24| ) . Such multi-phase similarity reductions are considered in chapter 4 of 
p6|, and we will discuss how our approach extends to these solutions in 



21 1 . Here we consider dependence on x and t n only, and henceforth we drop 
the suffix n, bearing in mind that the actual form of H depends on which 
particular flow we have chosen. 

Substituting the similarity forms into the equations of motion ( |2.1| ) and 
(|2.2| ) (and cancelling out powers of 9 on either side) yields the ODEs for w 
and y. If we let H denote the rescaled H (in terms of w with powers of 6 
divided out) then we obtain the following: 

(<9 3 + 8awd + Aaw') \ S w H[w] - ^-z^j = 0, (2.3) 
d(S y H[M[y]]+zy) = 0. (2.4) 
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The symbol d denotes derivatives with respect to z. Both equations ( |2.3| , [2.4|) 
can be integrated once, to yield 



d 2 f A 2 — (^-) 2 

_J_ + 4 awf + = 0, (2.5) 

$L-4ayf + \ = 0. (2.6) 

In the above, A is a constant of integration, and we have introduced the 
quantity 

/ := 5 W H - -U. (2.7) 



The equations ( |2.3| , |2.4|) deserve some comment, as we shall make use of 



them in section 3. With the definition ( |2.7| ), we may take / = f[w], and thus 
Q2.5Q is to be considered as the ODE for w, obtained directly from ( |2.3|) after 
an integration. However, using the scaled Miura map it is also consistent 
to take / = and so (by a slight abuse of notation) the equation 



(|2.6j ) represents the ODE for the variable y, obtained by integrating (|2.4fl 
However, a more obvious direct integration of (|2.4j ) would be 



5 y H + zy + a = 0, (21 



for some constant a. To show that ( |2.6| ) is equivalent to (|2.8|), it is necessary 
to use the scaled version of the identity 

5 V H = (M')*5 U H = (d x - Aav)5 u H 

(with M' being the Frechet derivative of M), and then it is found that the 
constants of integration are related by 

a = A . 

4a 

Another important point is that the same constant A appears in both (|2.5|) 
and Q2.6p . This is because the Miura map for the original PDEs ( |2.1| , ^.2[ ) be- 
comes a one-one correspondence between the ODEs for their scaling similarity 
solutions. The scaled Miura map is 

w = -y'- 2ay 2 , (2.9) 

and it has an inverse given by 

» = (2.10) 
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In ( |2.10| ) / should be regarded as f[w), and it is necessary to assume / ^ 
since otherwise this equation breaks down. It is straightforward to see the 
one-one correspondence between (|2.5|) and (|2.6|). Indeed, given a solution w 



of ( |2.5| ), the modified variable y may be defined by the inverse Miura map 
( 2.10| ). A direct calculation then shows that the Miura formula fl2.9|) holds, 



meaning that / can be reinterpreted as /[M[y]]. Thus (on rearranging ( p.lOp ) 
y satisfies the ODE (|2.6| ); the converse follows by reversing this argument (or 
by scaling the Miura map for the PDEs). To make things more concrete, it 
is worth looking at some particular cases. 

2.2 Similarity Solutions of KdV and PII 

The correspondence between PII (the equation (|1.3Q ) and scaling similarity 

h 

2 



solutions of the ordinary KdV equation is well-known. Taking H = ^u 2 in 



d3) with a = 1/2 yields KdV, 

u t = u xxx + Quu XJ (2.11) 
with scaling similarity solutions 

u(x,t) = (St)~*w(z), z = (St)-^x. 
The scaled Hamiltonian H = ^w 2 gives 

/M =w_ f> 

whence (|2.5| ) with a = | becomes 

„ o a(a + 1) + w' — (w') 2 

w" + 2w 2 -zw + — ^^ = 0, 2.12 

2w — z 

taking A = a + |. 

For the scaling similarity solutions of mKdV ( |1.2| ), the one-one correspon- 
dence between solutions of PII and ( p. 12 ) is given by 



w = -y' - y 2 

(the scaled Miura map) and 

w' + a 

y = o ' 

2w — z 

the latter being a particular case of the inverse Miura formula ([2.1 0|) . Observe 
that, in terms of y, we have 

f[M[y]] = -y'-y 2 -^ 

and on putting this into (|2.6|) with a = ~, PII results. This example is also 
considered in H, for instance. 
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2.3 Similarity Solutions of Fifth Order Equations 

Our second example constitutes the main subject of sections 3 and 4. Fol- 
lowing JTBf, we use the equation (|2.1|) with Hamiltonian density 

H = --ul - -bu 3 , 
2 x 3 

leading to the fifth order equation 

20 

ih = {u xxxx + (8a - 2b)uu xx - 2(a + b)u 2 x - —abu 3 ) x . (2.13) 
There are only three cases (i),(ii),(iii), corresponding to 

«4 i, j, »=~, -3, -4, (2.14) 

for which an equation of the form ( |2.13| ) is integrable, and these three cases 
correspond to the SK equation (|1.5| ), the fifth order KdV equation (|1.6| ), and 



the KK equation (|1 . 7|) respectively. When specializing to one of the three 
hierarchies, it is necessary to take the correct values of a and b in each case. 
While some of the calculations we present are valid for arbitrary a, b, all 
properties relevant to the integrability of the equations are lost in general. 
The scaling similarity solutions of fl2.13|) take the form 
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u(x, t) — (— 5t) 5 w(z), z = {—hi) 5£. 
From the scaled density H we find 

f\w] = w" — bw 2 — — , 
1 J 4a 

and substituting this into (|2.5|) yields a fourth order ODE for w, which we 
prefer not to present explicitly. In the next section we show that this ODE 
for w is equivalent to a non-autonomous Hamiltonian system of Henon-Heiles 
type. 

Scaling similarity solutions of the associated modified equations take the 
form 

v(x,t) = (-5t)~5y( Z ). 

Using the scaled Miura map, we may express / in terms of the modified 
variable y: 

f[M[y]] = -y'" - Aayy" - (4a + b){y'f - Aaby 2 y' - Aa 2 by A - f . (2.15) 

4a 
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Then the ODE for y is (from (gg)) 

= _ 2 ( 6a + 5) y y + 4a ( 4a _ 5) ( y y + y ^/) 2 ) + 16a 3 &2/ 5 + 2y + a? ( 2 . 16) 

with 

a = A . 

4a 

Given a solution to ([2.16D we can then obtain a solution to the fourth order 
ODE for w, via w = —y' — 2ay 2 (and vice- versa, using the inverse Miura map 



On substituting in the relevant values of a and b for the integrable cases 
(i),(ii),(iii) into the equation ( |2.16| ) we find 



y ™ = -5y'y" + 5(y 2 y" + y(y') 2 )-y 5 + zy + a, (2.17) 
yW = 10(y 2 y" + y(y') 2 )-6y 5 + zy + a, (2.18) 
yW = 5y'y" + 5(y 2 y" + y(y') 2 )-y 5 + zy + a, (2.19) 

the ODEs for scaling similarity solutions of the modified SK, fifth order 
mKdV and modified KK equations respectively. Notice that ( [2.17|) and (|2.19|) 



differ only by a sign in the even (y'y") terms. Hence if y^) is as solution to 
( |2.16j ) for case (i), then y^ = —y(i) will be a solution to that equation for 
case (iii) with a replaced by —a. This is because the modified hierarchies 
in these two cases are essentially the same; both SK and KK have a third 
order Lax operator, which factorizes to yield the Miura map (see [|15|]). In 
the next section we shall use this connection between case (i) and case (iii) 
to derive the Backhand transformations for these scaling similarity equations, 
after relating them to Henon-Heiles systems. 



Remark. In |24j the sequence of equations (|2.6| ) are referred to as the 
second Painleve equations of 2nth order (corresponding to scaling similar- 
ity solutions to the flow given by H n+ i). For instance, (up to a rescaling) 
the equation fl2.18j ) is the second Painleve equation of fourth order in the 
nomenclature of pi. However, this sequence of ODEs is by no means new; 



it appears in [14| under the name of the Painleve II Family. 



3 Henon-Heiles Systems 

3.1 Stationary Flows and Integrable Henon-Heiles 

The original Henon-Heiles system is given by a Hamiltonian with two degrees 
of freedom: 
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h = \{p\ + vl) + aqiql - -bqf, (3.1) 
The equations of motion are just Hamilton's equations, 

dq^ = dh_ dpj_ = _dh_ 

dz dpj ' dz dqj 

(We are denoting the time by z here to make connection with our other 
results.) It has been known for some time from Painleve analysis |H| that 



this system is integrable for only three values of the ratio r = a/b, namely 

r = -l, -1/6, -1/16. 
More recently it was shown by Fordy ]16| that in these integrable cases the 



equations of motion are just disguised versions of the stationary flows of 
the SK, fifth order KdV and KK equations. Thus the special values of a 
and b listed in the previous section give the right values for the ratio r in 
the cases (i),(ii),(iii) respectively. The zero curvature form of these PDEs 
yields a matrix Lax representation of the stationary flows, and then traces 
of powers of the Lax matrix give the Hamiltonian and a second constant of 
motion in involution (implying Liouville integrability) . It was subsequently 
shown that all three Hamiltonian systems are completely separable in suitable 
coordinates, and may be integrated in terms of theta functions of genus one 
(cases (i) & (iii)) or genus two (case (ii))||. 

3.2 Non- autonomous Hamiltonians for Scaling Simi- 
larity Solutions 

In order to consider the ODEs for scaling similarity solutions of the three in- 
tegrable fifth order PDEs (all of the form ( |2.13|) ), we adapt Fordy's approach 



161 to the stationary flows, and write them in Hamiltonian form. Recall that 



the previous section we have the quantity / given by 



f = w"- bw 2 - f , (3.3) 
4a 



and on substituting this into Q2-5|) a full fourth order ODE for w results. 
Instead of presenting this explicitly, we set 

w = qi, f = -aql, 



9 



the latter is a squared eigenfunction substitution [Q. Then we may rewrite 
(|2~5| , [3.3|) as a coupled system for g 1; q 2 : 



q" = bqi-aql + ^ (3.4) 

A 2 
4a^g2 



This coupled system just follows from Hamilton's equations ( |Q[ ), where now 
the Hamiltonian is given by 

1 1 A 2 1 

h = 2 (Pi + P2) + a 9i<?2 - 3^1 - ^2^ 2 ~ J-^i- (3-6) 



Compared with Q3.1|) , this has an extra inverse square term and a non- 
autonomous (time-dependent) term in the potential. 

Thus we have shown that the equations for the similarity solutions of 
(|2.13j ) may be viewed as non-autonomous Henon-Heiles systems. Because of 



the explicit time-dependence, the Hamiltonian h(z) satisfies a fourth order 
ODE which (in the generic situation /' ^ 0) is equivalent to the system 
(|3.4j,|3.5|) (this is proved in p(J). Applying Painleve analysis to this system 



(which is almost identical to the analysis for the autonomous case [10]), it 
is found that the Painleve test is only satisfied for the three special values 
of the ratio r (which correspond to the scaling similarity solutions of the 
integrable PDEs). Because they are similarity reductions of PDEs solvable 
by inverse scattering, these special cases can also be derived from a zero 



curvature condition p0| , and solved by the inverse monodromy method [14 
but we will not consider this approach here. 

For the particular special choices of a and b given in fl2.14j ), we denote the 
Hamiltonian systems (with Hamiltonian (|3.6|)) by H^), H(u), H(m) (corre- 
sponding to the reductions of SK, KdV and KK respectively). Each system 
has Backlund transformations which can be viewed as a canonical transfor- 
mations in the variables qj,pj,z, i.e. under the transformation 

(Qj,Pj,z) — * (Qj,Pj,z), 

the two-form 

uj = dpj A dqj — dh Adz (3.7) 

is preserved. 

The construction of the Backlund transformations is possible due to the 
existence of the Miura maps. For the system Titu) we have (by the general 
arguments of section 2) a one-one correspondence between its solutions and 
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the solutions of the modified similarity equation ( j2.18[ ), while both H(i) and 



TCuu\ have similar correspondences with Q2.17|) . We are also able to define a 
tau-function for the Hamiltonians, via 

h(Qj(z),Pj(z),z) oc (\og[T(z)}y 

(for some constant of proportionality, dependent on which of the three cases 
is being considered). We consider the system Tt(a) first, as this is the simplest, 
while the systems TC(i)jui) are then considered in tandem (because they are 
connected via the correspondence with the same modified equation ( [2.17|) ). 

3.3 The Hamiltonian System TLin\ 
The Hamiltonian for TC(u) is given explicitly by 

1 1 A 2 1 

hx = ~(p\ + vl) + ^ + 9i ~ y<? 2 ~ 2 " 2*0i- ( 3 - 8 ) 

We will henceforth put suffixes on all quantities to denote their dependence 
on A, as the Backlund transformation relates the same quantities for different 
values of this parameter. There is also the alternative parameter a = A — \ 
appearing in ( |2.18| ), but A is the more natural one in that ( |3.8| ) depends only 



on A 2 , i.e. 

h\ = h_ x , 

and thus the same is true for the solutions (qj(z),Pj(z)) to the system with 
Hamiltonian h\. The most convenient variables to use for this system are 

w\ = q u fx = --q 2 , 

where (by (|3.3|) with a = 1/2, b = —3) we have 

f x = w'{ + 3wl- Z -. (3.9) 

Using the Miura map, there is a one-one correspondence with solutions 
to ( p. 18 ). In other words, given a solution to the Hamiltonian system, we 



find a solution y — y\ to equation (|2.18|) (with a — A — |) via the formula 

CT), 

y> = 4^- ( 3 - 10 ) 

We also have the usual Miura expression, 

™x = -y'x-vl (3-ii) 
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which means that conversely a solution to Q2.18|) determines a solution to the 



system Hua. Before presenting the Backlund transformation, it is helpful to 
introduce the tau-function. 

Definition 1 For a solution (qj(z),pj(z)) of the system Tira), we have the 
Hamiltonian h\(z) = h\(qj(z),pj(z), z). The tau-function associated with 
this solution is given by 

h\{z) = -— log[r A (»]. 

The above definition is chosen to be consistent with the tau-function of 
the KdV hierarchy, where the dependent variable u is expressed as 

u{x,t) =2{log[T{x,t)}) xx . (3.12) 

For scaling similarity solutions we require that r depends on x and t only 
through the combination z = x6(t), which means we must have 

w x = 2(log[r(z)])". 
This agrees with our definition, for differentiating the Hamiltonian gives 

-(log[r(z)]) ff = h' x = = -~w x . 



Proposition 1 The Backlund transformation for the equation ( 2.18J may 
be written in the form 

yx+i = -yx + - r - (3.13) 
Jx 

Moreover, this induces a canonical transformation from the system H.(u) with 
parameter A to the same system with parameter A + 1: 

Qj -> Qj, Pj -> Pj, z ->• z, h x ^ hx+i- 

The modified variable y x may also be written in terms of the two tau-functions 
related by this transformation: 

y x = (log[T X ^/r x }y. (3.14) 

Proof. The first thing to observe is that w x is related to two different 
modified variables by the Miura map (|3.11|) : 



w x = -y' x -y\ = -yi x - y\. 
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This is a well-known property of Miura maps, but can be seen directly from 
the fact that the fourth order ODE for w\ just depends on A 2 , so 



W X = 10-A- 

The inverse Miura map (|3.10|) gives 



V±x = , (3.15) 



the solutions to ( |2.18| ) for parameter a = ±A — |. Looking at the modified 



equation (|2.18| ), we see that it is unchanged on sending 

y-*-y, « 

Hence — y_A will be a solution to this equation for a = A + ~ = (A + 1) — ~, 
or in other words we may identify 

2/A+l = -2/-A- 



It is then straightforward to derive ( |3.13|) using the inverse Miura formula 



This same argument also works for PII and the rest of the PII Family. Note 
that f\ in the right hand side of the Backlund transformation may be written 
in terms of y\ and its derivatives. 

To define the induced transformation of the Hamiltonian system H,(u), 
first of all it is necessary to write yx+i in terms of the variables appearing in 
hx- 

P2?2 + A 

2/A+l = 2 • ( 3 - 16 ) 

55 

The Miura map produces a new solution to the fourth order equation Q2.5|) 
at parameter value A + 1, 

^a+i = -y'x+x -yl+v 

Now we can define the new variables (<jj,pj) via 

qi = W\+i, pi=w'x +l , h = (-2/ A+ i)a, p 2 = -(-2/a+i) _5 /a+i> 



and these will clearly satisfy the system Hiu) with Hamiltonian hx+i- It is 
simple to write the new variables in terms of the old (and vice- versa), and 
we present the formulae here for completeness: 

9i = -gi-2T/A+i> 9i = ~f 1-2 ^A+i^ 

Pi = -Pi-4y\+i(.Qx + V\+i), Pi = 7|i + %a+i(9i + 2/a+i)> (317) 

q 2 = Ts, q 2 = fs, 



P2 = 2/A+1T2 + (A + 1)T 2, p 2 = - yx+1 r 2 - \r~ 
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where 



T = -ql + 8yx + ipi+8yi +iqi -4ql+2z, T = -q 2 2 -8y x+ m+8yl +1 qi -Aq{+2z 



and in the left-hand formulae yx+i is to be interpreted as the function of P2 
and q2 given by ( |3.16| ), while on the right-hand side take 



y\+i 



P2Q2 - (A + 1) 

?1 



For the sake of clarity, we present the (invertible) transformations in the 
variables w\ and y x diagramatically: 



yx+i 



(3.11) 



(3.11) 



( 13.171 ) 

W\ > Wx+l 



It is obvious that ( |3.17| ) is a canonical transformation, since the equations 
for both sets of variables are Hamiltonian; it is also possible to verify explicitly 
that the two-form fl3.7| ) is preserved. A direct calculation also shows that 
difference in the Hamiltonians is 



hx+i - hx = yx+x, 
which yields ( |3.14|) immediately. This completes the proof. □ 



(3.18) 



3.4 The Hamiltonian Systems TLu\ and TLim) 
The Hamiltonian for T~Lu) is 

1 1 1 \ 2 I 

U L /■ 2 , 2\ , L 2 , x 3 A -2 1 

x = - (Pi + P 2 ) + 2 qiq ^ + Q qi ' Y ?2 ~ 2 ZQl) 
while for 7Y(m) it is 

Hx 



(3.19) 



(3.20) 



To avoid confusion between the two, we use lower/upper case letters for the 
variables of the systems H/^/Tian) respectively. It is convenient to use the 
variables 



w x = qi, 



fx 



1 



W x = Qi, 



1 
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Hence, using ( |3.3| ) with the appropriate values of a and b, 



h = < + - f , ^ = w; + awI - z. 

For case (i) the Miura map and its inverse are given by 

= -y' x - vl yx = ^t^. ( 3 - 21 : 



where y\ is a solution to ( |2.17| ) for a = A — |. Although case (iii) is related 



to the same modified equation, it will be helpful for deriving the Backlund 
transformation to return to the original formalism of Section 2, where the 
case (iii) (with a = 1/4) the Miura map and its inverse are given by 

W X = -Y{ - \Yl Y X = (3.22) 



where Y\ is a solution to ( |2.19| ) for a = A — 1. The derivation of the Backlund 



transformation for the equation (|2.17|) (or equivalently (|2.19 )) is most easily 



achieved with the tau-functions, naturally related to the tau-functions of the 
SK/KK hierarchies (see |2l| for definitions of these). 



Definition 2 For a solution (qj(z),pj(z)) of the system Ti^), with the Hamil- 
tonian h\(z) = h\(qj(z),pj(z), z) , the associated tau-function is given by 

M*) = -3^1og[T A (*)]. 



Definition 3 For a solution (Qj(z), Pj{z)) of the system 'Hun), with the 
Hamiltonian H\(z) = H\(Qj(z), Pj(z), z), the associated tau-function is given 
by 

^) = -~log[7>(z)]. 

Now we may show: 



Proposition 2 The Backlund transformation for the equation (\2. 1 Tj) may 

be written in the form 



A 2(A + 

yx+3 = y\- -j- h — 



(3.23) 
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This is related to a canonical transformation between the system Ti^ with 
parameter A and the system Ti(iu) with parameter A — | : 

qj^Qj, pj -> Pj, z -> z, h x ^>H x _s. 

The modified variable y\ may also be written in terms of the two tau-functions 
related by this transformation: 

y x = (log[f A _3/r A 2 ])'. (3.24) 

Proof. As for case (ii), w x (and also W x ) may be related to two different 
modified variables by the Miura map. Hence we see that 

yx - y-x = 7-, (3.25) 
Jx 

2A 

Y x -Y_ x = — . (3.26) 
r\ 

Clearly ( |3.25|) constitutes a Backlund transformation for ( |2.17 ), as it relates 
two solutions for different parameter values. However, it is not very useful 
because it cannot be iterated to obtain a sequence of solutions. To achieve 
this requires a canonical tranformation from 7iu) to Tirm), and then from 
TL(iu) back to Ti(€\ with the overall change A — > A + 3. First of all observe 
that on comparing the parameters a in ( |2.17|) and (|2.19|) , it is apparent that 
we may identify 

yx = -Y-\+§, 

and so the Miura map for case (iii) implies 

W„ x+ 3 =y' x - -yl 
Since W_ x+ 3 may be found from w x , via the formula 

W_ x+I = -w x - *~yl (3.27) 

it is obvious that there is an induced canonical transformation from the 
system with Hamiltonian h x to the system with Hamiltonian H_ x+ a = H x _3. 
This transformation and its inverse may be calculated explicitly in terms of 
coordinates, i.e. 



Q\ = -qi ~ \vl, qi = ~Qi - Ivl, 



Pi = -Pi + 3y x {qi + y 2 x ), pi = -Pi-tyxiQi + ivl), 
Q 2 = Ks, q 2 = & z 

P2 = -hx^ -2(\-m- 1 2, p 2 = y x & + \X 



~ i 

2 . 
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where 

N = -2q 2 2 - \2y xVl - Qq\ + 6z, N = ~Q\ + 6y x Pi - 3Qj + 9ylQi + -^y\ + 3z, 

and in the above y\ is to be interpreted alternately as a function of p 2 , q 2 
and P 2 , Q 2 : 



tjx 



p 2 q 2 - A 2P 2 Q 2 + 4(A 



4 



Ql 



This is the analogue of the transformation between autonomous Henon-Heiles 
systems considered in [|, |I~3|j . 

The tau-function formula fl3.24j ) follows directly from a calculation of the 
difference of the two Hamiltonians: 



h x - H x _z = -y x . 



Then in terms of tau-f unctions, we have 

d 



.'/a+:; - .'/,\ = j~ z l lo g[ f A+|/^A-|] + 21og[r A /r A+3 ] 



On using the formulae ( |3.25| , |3~2"(3D , the Backlund transformation fl3.23| ) fol- 
lows. Thus overall there is an induced canonical transformation from Ti^ to 
itself, with the parameter A — > A + 3, as can be seen from a diagram: 



(3.21) 



y-x 



W-X 



-Y x+ I > Y-x-l 



-yx+3 



(3.27) 



MA 



^A+| 



( P-22| ) 
: W 



(3.27) 



( |3.2l|) 



Note that in the right hand side of ( |3.23|) , f\ may be determined entirely in 
terms of yx and its derivatives, and similarly for F x+ 3. □ 



3.5 Analogues of the Toda Lattice for Sequences of 
Tau-Funct ions 



It is well known p3| , p7| that for PII, the tau-functions related by the Backlund 
transformation (after rescaling) satisfy the Toda lattice equation, 



Dlrx ■ T X 



Tx-lTx+l- 



(3.28) 



The tau-functions of Q2.18] ) (related by the Backlund transformation) also 
satisfy a bilinear lattice equation, which is a fourth order analogue of G3.28| ). 
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Indeed, on the one hand (substituting for w x in terms of r A into (|3.9|)) we 
have 



h = r- 2 (Dtr x .r x )- Z -, (3.29) 



while on the other hand (from (|3.14j ) and (|3.15|) ) 



(\og[r x+1 T X ^/r 2 x })' = y x + y_ A = (log[/ A ])'. (3.30) 



Thus integrating left- and right-hand sides of ( |3.30| ) and comparing with 
(|3.29| ) yields the bilinear equation 



2D z t x ■ t x - zr x = kxTx-iTx+v (3.31) 

Note that the constant of integration k x may be rescaled arbitrarily, since 
the tau-functions can always be rescaled (without affecting h x ,y x ). We take 
the convention k x = 1, which ensures that the tau-functions for the rational 
solutions (presented in section 4) are monic polynomials in z. However, 
the derivation of ( |3.31| ) made the generic assumption f x ^ 0, This may be 
violated for A = 0. In this case the constant in ( p.31|) vanishes (i.e. ko = 0), 
and solutions related to PI (discussed in the next section) are obtained. 

Almost identical arguments lead to the following equations for the tau- 
functions of the systems Ti^ and 7i(iuy. 

™i - 6t 2 x D 4 z t x ■ r x + 9(D 2 z r x ■ r A ) 2 = f A _|f A+f , (3.32) 

3 

~ 1 D % ■ f A = r A _|r A+ |f A . (3.33) 

We have used the same conventions and genericity assumptions as for H,(u). 
Observe that these equations are no longer bilinear, but it is consistent to 
assign weight one to r and weight two to f in the two equations. 



4 Special Solutions 
4.1 Rational Solutions 

The rational solutions of the Painleve equations have been constructed in 
many different ways || g, |23[. The rational solutions of Q2.18D follow a 



pattern very similar to those of PII. They may be generated by applying 
the Backlund transformation ( |3.13|) repeatedly, starting with the solution 



yi = 0, to produce a sequence of solutions y x for half-integer values of A. It 
is obvious that the solutions thus generated must be rational, as the Backlund 
transformation only involves differentiations and algebraic operations. The 
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Miura map (|3.11| ) also yields the associated sequence of rational functions 
w\ (giving solutions to the system TCua). From an algorithmic point of view, 
it is much more convenient to calculate the rational solutions from their 
polynomial tau-functions T\, which can be obtained by iteratively solving 
( |3.31| ). We present a few of these rational solutions, and their associated tau- 



functions, in table 1. Note that it is sufficient to consider only the positive 
half-integer values of A, because of the formula y\ + i = —y~\. 



A 


1/2 


3/2 


5/2 


7/2 


9/2 


y\ 





_ i 


_2 


-3(z 5 +96) 


-4(z 15 -72z 10 +217728z 5 -1741824) 


z 


z 


z(z 5 -144) 


z(z 15 -1152(z 10 -84z 5 +6048)) 


W\ 





2 


6 


-12(z 10 +432(z 5 +8)) 


-20z 3 (z 15 +1008z 10 +943488z 5 -47542144) 


z 2 


z 2 


z 2 (z 5 -144) 2 


(z 10 -1008(z 5 +48)) 2 


T\ 


1 


Z 


Z 3 


z 6 - UAz 


z w - 1008(^ 5 + 48) 



Table 1: Rational solutions for case (if). 



The solutions above can also be obtained as specializations of the well- 
known sequence of rational solutions of the mKdV/KdV hierarchy constructed 
in f|], and thus it is possible to express the polynomial tau-functions r A as 
Wronskian determinants [|2D| . 

Rational solutions y\ of the equation ( |2.17| ), together with the Miura- 
related w\, W x _3 (giving solutions to systems H^^m)) may be generated in 



a similar fashion, by applying the Backlund transformation (|3.23|) . Some of 
these solutions are given in table 2. There is one sequence for the values 
A = 3j + \ (for integer j), generated from the trivial solution yi = 0, and 

a second sequence for A = 3j — \ (which is related to the first by (|3.25| )). 
Once again, it is easier to generate solutions using the equations ( |3.32| , |3.33| ) 
for the tau-functions. We present some of these in table 3 below (using two 
rows for the pair f A± 3 related by ( |3.32| )). 



4.2 Solutions Related to PI 

All three systems have special solutions related to ( [L.l|) , the first Painleve 
transcendent (PI). If we consider the system Ti^ in the case A = 0, the 
same substitution that works in the ordinary (autonomous) system causes 
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A 


-7/2 


-5/2 


-1/2 


1/2 


5/2 


7/2 


11/2 


V\ 


4 


3(2 5 -24) 


i 





_2 


_3 


-5z 4 (2 5 +216) 


2 


2(2 5 +36) 


2 


2 


z 


(2 5 +36)(2 5 -144) 


W\ 


12 

2 2 


6 

2 2 








6 

2 2 


12 

2 2 


-302 3 (2 5 +576) 
(2 5 -144) 2 


W\ 3 


12 

2 2 


-152 3 (2 5 -144) 
2(2 5 +36) 2 


3 

22 2 








3 
2z 2 


12 

2 2 



Table 2: Rational solutions for cases (i) and (iii). 



A 


1/2 


5/2 


7/2 


11/2 


T\ 


1 


z 


z 2 


^ 5 - 144 




1 


1 


z 


z 5 + 36 




z 


z 5 + 36 


z s 


z ie - 2 7 .3V 1 + 2 U .3 4 .11^ 6 + 2 14 .3 6 .llz 



Table 3: Polynomial tau-functions for cases (i) and (iii). 

the equations of motion to separate. Putting 

Q± = qi± <? 2 
into Hamilton's equations for ho, we find 

<y± = -\?±+\, 

which (up to a scaling) is just two separate copies of PI. The corresponding 
solution to ( |2.17| ) is 

y = (log[Q + - Q-})', 

where we assume that Q + and Q_ are not equal. So by applying the Backlund 
transformation ( |3.23| ) we get the general solution to the system Ti^ for A = 
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3j, and to H(iu) for A = 3(j + |), for all integers j. However, there is also 
the degenerate case Q + = for which 

/o = 0. 

This implies that the inverse Miura map and the Backlund transformation 
both break down. However, it is still possible to obtain a sequence of special 
solutions for the same parameter values, and they are also related to PI. 
Similarly, starting from a degenerate solution Y , corresponding to 







the Backlund transformation ( p.23|) gives a sequence of special solutions to 



the system Ti^ for A = 3(j + |), and to H(m) for A = 3j. We explain these 
degenerate solutions in more detail for case (ii). 

The degenerate case for all three systems is A = 0, f = 0. So in case (ii) 
w = w must satisfy 



w" + 3w 2 -- = 0, 

which is equivalent to PI (after rescaling w and z). The inverse Miura map 
breaks down in this case. However, the ordinary Miura map means that 
y = yo satisfies the Riccati equation 

y' + y 2 + w = 0. 

This is linearized by setting y = (log[r])', to yield 

t" + W T = 0, 

which (for each wq) gives a one parameter family of solutions yo to fl2.18| ). 



z 



Because fo = 0, the Backlund transformation ( 3.13 ) breaks down for the 
solution y . However, we still have 

yx+i = -y-x, 

and so we can safely apply the Backlund transformation to yi = —yo to 
obtain a sequence of solutions for all integers A, as well as corresponding 
solutions to the system TC(u). 

Remark. These degenerate solutions to (|2.18|) are the analogues of the 



Airy function solutions to PII (see [23] and some of the other references). 



The fact that solutions of PI give special solutions of ( |2.18| ) is a particular 



case of a result in 24 . 
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5 Conclusions 



We have found that the Hamiltonian formalism for stationary flows of cer- 
tain integrable PDEs may be extended to the similarity solutions. For the 
particular examples we have considered, non-autonomous generalizations of 
integrable Henon-Heiles systems are obtained. Applying the same technique 
to the travelling wave similarity solutions of ( |2.13| ) leads to Henon-Heiles 
Hamiltonians with harmonic terms, while the scaling similarity solutions of 
Hirota-Satsuma and a related system correspond to non-autonomous versions 
of the integrable Hamiltonians with quartic potentials considered in ]7], §; 
these examples appear in PU| . Our method should always work provided 
that the ODEs for the similarity solutions can be written in the form 



B w f = 



where B w is the scaled version of the Hamiltonian operator for the PDE, and 
/ is the gradient of the (scaled) Hamiltonian for the PDE plus some extra 
(possibly non-autonomous) terms (cf. section 2). 

It seems to be a common phenomenon that equations of Painleve type 
can be viewed as non-autonomous generalizations of autonomous, integrable 



Hamiltonian systems. For instance, the equations in [jn| are non-autonomous 
versions of the integrable Neumann systems. This connection might be used 
to study the asymptotics of higher order Painleve equations. Recently we 
have succeeded in extending our approach to the multi-phase similarity so- 



lutions of the kind considered by Newell [26|, but this will be considered 



(together with the connection to isomonodromic deformations) in a forth- 



coming article 21 . 
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